
function y=funcs(theta)
global tauk taun  alpha fgdp phi b delta Time
%equation that matches the formal gdp
y=fgdp-theta*((((((1-tauk)*theta*alpha)/(1/b-1+delta))^(1/(1-alpha)))*(Time/(1+phi+((alpha*(1-tauk)*phi)/((1-alpha)*(1-taun)))*((1-b)/(1-b+delta*b)))))^alpha)*((Time/(1+phi+((alpha*(1-tauk)*phi)/((1-alpha)*(1-taun)))*((1-b)/(1-b+delta*b))))^(1-alpha));
end